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ABSTRACT 

The combined analysis of the cosmic microwave background on large scales 
and Lyman-a forest on small scales provides a sufficiently long lever arm to 
obtain strong constraints on the slope and curvature of the power spectrum 
of primordial density fluctuations. We present results from the combination of 
the first year WMAP data and the dark matter power spectrum inferred by 
Viel et al. (2004) for two different sets of high resolution and high signal-to- 
noise quasar absorption spectra: the Croft et al. (2002) sample with a median 
redshift z = 2.72 and the LUQAS sample (Kim et al. 2004) with a median 
redshift z — 2.125. The best fit value for the rms fluctuation amplitude of 
matter fluctuations is as = 0.94±0.08 and n s = 0.99±0.03, if we do not include 
running of the spectral index. The best fit model with a running spectral index 
has parameters n s = 0.959 ±0.036 and n run = —0.033 ±0.025. The data is thus 
consistent with a scale-free primordial power spectrum with no running of the 
spectral index. We further include tensor modes and constrain the slow-roll 
parameters of inflation. 

Key words: Cosmology: observations - cosmology: theory - cosmic mi- 
crowave background - quasars: absorption lines 



1 INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe team 
(WMAP, Bennet et al. 2003, Spergel et al. 2003) has 
presented impressive confirmation for what is now con- 
sidered the standard cosmological model: a flat uni- 
verse composed of matter, baryons and vacuum en- 
ergy with a nearly scale-invariant spectrum of primor- 
dial fluctuations. One of the surprising results of the 
WMAP team's combined analysis of cosmic microwave 
background (CMB), galaxy redshift survey and Lyman- 
a forest data was the evidence for a curvature of the 
power spectrum of primordial density fluctuations. This 
was expressed by the WMAP team in terms of a non- 
vanishing derivative n run = dn s /d\nk of the spectral 
index n s . In inflationary models the shape of the primor- 
dial power spectrum is directly related to the potential 
of the inflaton field. For slow-roll inflation simple re- 
lations exist between the slow-roll parameters and the 
ratio of tensor to scalar fluctuation, the power-law index 
and the derivative of the power-law index of primordial 
fluctuations (see Liddle and Lyth 2000 for a review). 
Constraints on the power spectrum of primordial den- 
sity fluctuations can thus at least in principle constrain 
inflationary models (Hannestad et al. 2002, Peiris et al. 



2003, Leach & Liddle 2003, Kinney et al. 2004, Tegmark 
et al. 2004). 

The Lyman-a forest data constrains the dark mat- 
ter (DM) power spectrum on scales of ~ 4 — 40 Mpc//i 
in the redshift range 2 < z < 4 (Croft et al. 1998, Hui 
1999, Croft et al. 1999, McDonald et al. 2000, Hui et al. 
2001, Croft et al. 2002 [C02], McDonald 2003, Viel et al. 
2003, Viel, Haehnelt & Springel 2004 [VHS]). Due to the 
long lever arm between Lyman-a forest on small scales 
and the CMB on large scales the combined analysis of 
these data sets alone can provide strong constraints on 
the slope and curvature of the primordial power spec- 
trum in the context of CDM-like models. Croft et al. 
(1999) inferred an amplitude and slope which was con- 
sistent with a COBE normalised ACDM model with a 
primordial scale invariant fluctuation spectrum (Phillips 
et al. 2001). McDonald et al. (2000) and C02, using a 
larger sample of better quality data, found a somewhat 
shallower slope and smaller fluctuation amplitude. This 
latter data was part of the WMAP team's analysis men- 
tioned above which gave evidence for a tilted primordial 
CMB-normalised fluctuation spectrum (n s < 1) and/or 
a running spectral index (Bennet et al. 2003; Spergel et 
al. 2003; Verde et al. 2003). There are, however, a num- 
ber of systematic uncertainties in the amplitude of the 
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DM power spectrum inferred form the Lyman-a forest 
data (C02, Zaldarriaga, Scoccimarro & Hui, 2003; Zal- 
darriaga, Hui & Tegmark, 2001; Gnedin & Hamilton 
2002; Seljak, McDonald & Makarov 2003, VHS). 

VHS have recently used the method developed by 
C02 together with a large suite of new high-resolution 
numerical hydro-simulations to obtain an improved es- 
timate of the DM power spectrum from the Lyman- 
a forest. Two large samples of high-resolution spectra 
were used for this analysis, the sample of HIRES/LRIS 
spectra in the original analysis of C02 and the Large 
Sample of UVES QSO Absorption Spectra (LUQAS) 
of Kim et al. (2004). VHS demonstrated that the in- 
ferred rms fluctuation amplitude of the matter density 
is about 20% higher than that inferred by C02 if the 
best possible estimate for the effective optical depth, 
as obtained from high-resolution absorption spectra, is 
used. VHS further found that with this effective optical 
depth the inferred DM spectrum is consistent with a 
scale invariant COBE normalized DM power spectrum 
thus confirming the findings of Seljak et al. (2003). The 
running spectral index model advocated by the WMAP 
team falls significantly below the DM spectrum inferred 
from the Lyman-a forest data (VHS). 

In this paper, we will perform a combined analysis 
of the results of VHS based on the two independent data 
sets of C02 and K04 together with the WMAP data. The 
plan of the paper is as follows. In Section 2 we present 
the two data sets and briefly outline the method used 
to infer the linear dark matter power spectrum. Section 
3 contains a description of the method used to estimate 
rms fluctuation amplitude as, the spectral index n s , a 
possible running of the spectral index and the param- 
eters of slow-roll inflationary models. In Section 4 we 
summarise and discuss our results. 



2 THE WMAP AND Lyman-a FOREST 
DATA 

2.1 The WMAP data 

The CMB power spectrum has been measured by the 
WMAP team over a large range of multipoles (I < 
800) to unprecedented precision on the full sky (Hin- 
shaw et al. 2003, Kogut et al. 2003). For the analy- 
sis we use the one year release of the WMAP tem- 
perature and temperature-polarization cross-correlation 
power spectrum 1 . The CMB power spectrum depends 
on the underlying primordial fluctuation via the pho- 
ton transfer functions. The angular power spectrum in 
the temperature exhibits acoustic peaks at I — 220 
with an amplitude of AT ~ 75 fiK and I = 540 with 
AT « 49 fiK (Page et al. 2003). The WMAP team 
(Spergel et al. 2003) showed that the measured CMB 
power spectrum is in excellent agreement with a flat 
universe and throughout this paper we will adopt this 
as a prior. In order to estimate parameters from the 
CMB we exploit for COSMOMC 2 an adapted version 

1 http://lambda.gfsc.nasa.gov 

2 http://www.cosmologist.info 



of the CMB likelihood code provided by the WMAP 
team (Verde et al. 2003, Hinshaw et al. 2003, Kogut et 
al. 2003). 

2.2 The Lyman-a forest data 

The power spectrum of the observed flux in high- 
resolution Lyman-a forest data constrains the DM 
power spectrum on scales of 0.003 s/km < k < 
0.03 s/km. At larger scales the errors due to uncertain- 
ties in fitting a continuum to the absorption spectra 
become large while at smaller scales the contribution of 
metal absorption systems becomes dominant (see e.g. 
Kim et al. (2004) and McDonald et al. (2004) for discus- 
sions). We will use here the DM power spectrum which 
VHS inferred for the C02 sample (Croft et al. 2002) 
and the LUQAS sample (Kim et al. 2004). The C02 
sample consist of 30 Keck HIRES spectra and 23 Keck 
LRIS spectra and has a median redshift of z = 2.72. 
The LUQAS sample contains 27 spectra taken with 
the UVES spectrograph and has a median redshift of 
z = 2.125. The resolution of the spectra is 6 km/s, 8 
km/s and 130 km/s for the UVES, HIRES and LRIS 
spectra, respectively. The S/N per resolution element is 
typically 30-50. 

The method used to infer the dark matter power 
spectrum from the observed flux has been proposed 
by Croft et al. (1998) and has been improved by C02, 
Gnedin & Hamilton (2002) and VHS. The method relies 
on hydro-dynamical simulations to calibrate a bias func- 
tion between flux and matter power spectrum: Pf(&0 = 
b 2 (k) P(k), on the range of scales of interest, and has 
been found to be robust (Gnedin & Hamilton 2002). 

There is a number of systematic uncertainties and 
statistical errors which affect the inferred power spec- 
trum and an extensive discussion can be found in C02, 
Gnedin & Hamilton (2002) ; Zaldarriaga, Scoccimarro & 
Hui (2003) and VHS. VHS estimated the uncertainty of 
the rms fluctuation amplitude of matter fluctuation to 
be 14.5 %. There is a wide range of factors contribut- 
ing to this uncertainty. In the following we present a 
brief summary. The effective optical depth which is used 
to calibrate numerical simulation and has to be deter- 
mined separately from the absorption spectra is a ma- 
jor systematic uncertainty. Unfortunately, as discussed 
in VHS there is a considerable spread in the measure- 
ment of the effective optical depth in the literature. De- 
terminations from low-resolution low S/N spectra give 
systematically higher values than high-resolution high 
S/N spectra. There is little doubt that the lower val- 
ues from high-resolution high S/N spectra are appro- 
priate and the range suggested in VHS leads to a 8% 
uncertainty in the rms fluctuation amplitude of the mat- 
ter density field (Table 5 in VHS). Other uncertainties 
are the slope and normalization of the temperature- 
density relation of the absorbing gas which is usually 
parametrised as T = To (1 + 5) 7_1 . To and 7 together 
contribute up to 5% to the error of the inferred fluctua- 
tion amplitude. Our estimate of the uncertainties due to 
the method was also 5%. Somewhat surprisingly VHS 
found substantial differences between numerical simula- 
tions of different authors and assigned a corresponding 
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uncertainty of 8%. Other possible systematic errors in- 
clude fluctuations of the UV background, temperature 
fluctuation in the IGM due to late helium reionization 
and the effect of galactic on the surrounding IGM. In 
VHS we had assigned an uncertainty of 5% for these. 
Summed in quadrature this led to the estimate of the 
overall uncertainty of 14.5% in the rms fluctuation am- 
plitude of the matter density field. 

For our analysis we take the inferred DM power 
spectrum in the range 0.003 s/km < k < 0.03 s/km as 
given in Table 4 of VHS. Note that we have reduced 
the values by 7% to mimick a model with 7 = 1.3, the 
middle of the plausible range for 7 (Ricotti et al. 2000, 
Schaye et al. 2000) . In the next section we will also show 
results for the DM power spectrum inferred by C02 in 
their original analysis (their Table 3). 



3 CONSTRAINING THE PRIMORDIAL 
POWER SPECTRUM AND INFLATION 

3.1 The Power Spectrum 

We will first attempt to constrain the primordial power 
spectrum of scalar fluctuations. We parametrise the 
power spectrum by 



n s (k) 



(1) 



where we choose the pivot scale to be fco = 0.05 Mpc -1 
for the analysis in this subsection. The simplest case is a 
scale invariant power spectrum with n s = 1 constant. In 
a more general context, and assuming d 2 n s /dlnfc 2 = 0, 
the spectral index is defined by: 



n 3 (k) = n s (ko) + 



dns 
dink 



In 



(fee) ' 



(2) 



with n run = dris/d lnfe parametrising the evolution of 
the spectral index to first order. We extended the COS- 
MOMC package (Lewis & Bridle 2002) to include the 
Lyman-a data as discussed in Section 2. As mentioned 
before the main uncertainty of the Lyman-a data is its 
relative calibration with respect to the underlying dark 
matter power spectrum. This is mainly influenced by the 
effective optical depth r e ff, the normalization and slope 
of the temperature-density relation, To and 7 and the 
other uncertainties described in section 2.2. We describe 
this uncertainty with an effective calibration amplitude 
A and obtain for the log-likelihood 



- In Chya (X ^ 
k 



[APdat(fc)- Pth(fe)] 2 



(3) 



We introduce the calibration A as an extra parameter 
into the Markov chain and multiply the likelihood with 
the prior 



Pprior OC exp 



1 (A-Ay 



(4) 



In order to obtain constraints on the cosmological pa- 
rameters we will marginalize over the calibration param- 
eter A. Note that this differs from the method described 



in Verde et al. (2003) . As a standard prior on the mean 
we choose A = 1. 

The cosmological parameters which are varied for 
the combined analysis are the physical baryon density 
fib h 2 , the physical dark matter density fi c h 2 , the an- 
gular diameter distance to last scattering 9, the optical 
depth r due to instantaneous complete reionization, the 
amplitude A s of primordial perturbations and as men- 
tioned above the spectral index of scalar perturbations, 
either running or constant. Note that we restrict the 
analysis to a flat universe with a cosmological constant, 
where we neglect the contribution from neutrinos and 
tensor modes. We further restrict the analysis to opti- 
cal depths r < 0.3 as suggested by the WMAP team 
(Spergel et al. 2003, Verde et al. 2003). The amplitude 
of the matter fluctuations in a 8 Mpc/h sphere, erg, is a 
derived quantity for this analysis. 

We first discuss the results when we assume a con- 
stant spectral index n s (n run = 0) . 

In Figure 1 we show the joint likelihood for the 
scalar spectral index n s and erg. On the left is the re- 
sult for WMAP only. Note that we only include the 
WMAP data set and not the extended CMB data set 
in order to concentrate on the constraints which are 
obtained by adding the Lyman-a forest. The WMAP 
data alone gives only very weak limits on the ampli- 
tude of the matter fluctuations and the spectral in- 
dex. For a marginalized likelihood the spectral index 
is n s = 0.990 ± 0.039, consistent with a scale invariant 
power spectrum as found by the WMAP team and other 
authors. In the middle panel we show the result if we 
combine the WMAP data with Lyman-a data from C02 
as compiled by Gnedin & Hamilton (2002) . Note that we 
impose an uncertainty on the calibration of a a = 0.25 if 
A < 1 and a A = 0.29 if A > 1 (Verde et al. 2003). The 
spectral index differs now from a scale invariant n s = 1 
primordial spectrum at the 2a level. The marginalized 
best fit value is n s — 0.955 ± 0.023. The right panel 
shows result for the WMAP data combined with the 
DM power spectrum as inferred by VHS. For this data 
we use a prior with a a = 0.29 throughout, which is 
twice the uncertainty of the rms fluctuation amplitude 
as discussed in Section 2. As discussed in detail by VHS, 
the inferred amplitude is about 20% larger compared to 
that inferred by the C02 data with the lower effective 
optical depth. The spectral index is consistent with a 
scale invariant spectrum with n s — 0.992 ±0.032. There 
is no significant evidence for a non-scale invariant spec- 
trum in agreement with the simpler analysis in VHS 
which obtained n s = 1.01 ± 0.02(stat.) ± 0.06(syst.) and 
cr 8 = 0.93 ± 0.03(stat.) ± 0.09(syst.). 

In Figure 3 we show the dependence of the as in- 
ferred from the combined analysis on the amplitude of 
the power spectrum inferred from the Lyman-a forest 
(parameterized as A, eq. 3). The dependence around the 
best fitting value is quite close to linearity (erf ~ A). 
This plot was obtained with a prior on A given by eq. 
(3). The contours close at small and large values due to 
the constraints by the WMAP data alone. On the same 
axis we show explicitly the dependence on the effec- 
tive optical depth (relative to the assumed best fitting 
value). The degeneracy of the inferred amplitude with 
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Figure 1. Joint 1- and 2-cr likelihoods on erg and n s . In the left panel the constraints from WMAP only, the middle panel from 
WMAP combined with the Lyman-a power spectrum of C02, for a T c fj = 0.349. On the right the constraints from WMAP 
combined with VHS, for a T e g(z = 2.125) = 0.17 and T ef j (z = 2.72) = 0.305. The dashed lines indicate the spectral index for a 
scale invariant n s = 1 primordial spectrum. 
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Figure 2. Joint 1- and 2-cr likelihoods on n run and n s . In the left panel the constraints are from the WMAP only, the middle 
panel from WMAP combined with the Lyman-a power spectrum of C02, for a r c g(z = 2.72) = 0.349. On the right the 
constraints from WMAP combined with the results of VHS, for a t cB (z = 2.125) = 0.17 and r cff (z = 2.72) = 0.305. The dashed 
lines indicate a scale-invariant n s = 1 and n rlm = primordial spectrum. 



the assumed effective optical depth already discussed by 
C02 and Seljak et al. (2003) is clearly seen. For the re- 
lation between amplitude A and effective optical depth 
we assumed A oc r~^ A as found by VHS. 

We will now discuss the constraints which we ob- 
tain if we include a running spectral index. In Figure 2 
we show the joint la and 2a likelihoods in the n run -n s 
plane marginalized over all other parameters. In the left 
panel we show the constraints from WMAP only. Again 
there is no significant evidence for a deviation from 
a simple scale invariant power spectrum, although the 
marginalized best fit values are n s = 0.922 ± 0.062 and 
n r un = —0.050 ±0.037. In the middle we show the result 
from combining WMAP with the C02 data. The best fit 
values are n s = 0.924T0.030 and n run = -0.052T0.016. 
This would rule out the most simple scale invariant 
power spectrum on the 3-cr level. This is consistent with 
the findings of the WMAP team but somewhat unusual 
for slow-roll inflationary models which require n run to 
be of the order of (|n s — 1|) 2 (Kosowsky & Turner 1995; 
but see Dodelson & Stuart 2002). If we instead com- 
bine the VHS data with WMAP (right panel) the con- 



tours shift again towards a scale invariant spectrum 
with no running of the spectral index. The marginal- 
ized one dimensional errors are n s = 0.959 ± 0.035 and 
^run — -0.033 ± 0.025. This makes the difference to a 
simple scale invariant primordial power spectrum sta- 
tistically insignificant. 

In Table 1 we show the constraints on all the cos- 
mological parameters for the combined WMAP and 
Lyman-a analysis with and without a running spectral 
index. 



3.2 Constraining Slow-Roll Inflation 

In order to obtain model independent constraints on in- 
flation we concentrate on slow roll models. Since we are 
using a running of the spectral index of the scalar per- 
turbations over a wide range in wavenumbers we use a 
second order extension of the standard slow-roll approx- 
imation (Liddle & Lyth 1993). This is described by the 
parameters (Martin & Schwarz 2000, Leach et al. 2002) 
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Table 1. The marginalized results on cosmological parameters. The quoted values correspond to the peaks of the marginalized 
probability distributions. Errors are the 68% confidence limits. 





WMAP+C02 (power law) 


WMAP+VHS (power law) 


WMAP+C02 (running) 


WMAP+VHS (running) 


n c h 2 


0.109 ±0.010 


0.124 ±0.013 


0.124 ±0.013 


0.132 ±0.015 


n h h 2 


0.022 ± 0.009 


0.024 ±0.001 


0.022 ±0.001 


0.023 ± 0.001 


h 


0.721 ± 0.040 


0.704 ± 0.047 


0.683 ± 0.048 


0.677 ± 0.047 


T 


0.115 ± 0.056 


0.157 ±0.068 


0.192 ±0.065 


0.182 ± 0.069 


<?8 


0.794 ± 0.052 


0.936 ± 0.079 


0.921 ±0.071 


0.978 ± 0.090 


n B 


0.955 ± 0.023 


0.992 ± 0.032 


0.924 ± 0.030 


0.959 ± 0.036 


TItuxi 






-0.052 ± 0.016 


-0.033 ±0.025 
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Figure 3. Joint 1- and 2-<x likelihoods on as as a function of 
the amplitude of the power spectrum A and of the effective 
optical depth. 
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where primes denote derivatives with respect to the 
field and M-pi = (87rG)~ 1//2 is the reduced Planck mass. 
Note that e\, €2, £3 correspond to the conventional slow- 
roll parameters by e\ = ev, £2 = 2ey — 2-qv and 
£2^3 = 4ey — Gevr)v + 2£v- If we calculate the spec- 
trum of primordial perturbations in the inflaton field 
we can identify the slow-roll parameters with param- 
eters of the power spectrum. In general inflation also 
predicts tensor fluctuations or gravitational waves. We 
neglect here the running of the tensor index which is 
given by nl un = e±e2- With the slow-roll conditions we 
obtain 



1 



-2ei — £2 — 2ex — Ce263 



(8) 
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Figure 4. Joint likelihood in the rjy-^y plane, where the 
colour coding of the dots is scaled with ey The contours are 
the 68% and 95% confidence levels. 
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(9) 
(10) 



with C « 0.73. We then obtain as a "consistency" rela- 
tion between the scalar to tensor ratio r and the tensor 
spectral index 



r r 



128 



(11) 



Note that this is only true if we neglect the running 
of the tensor spectral index. We will calculate the slow- 
roll quantities from the Markov chains generated in the 
parameters n s , n run and r. In this case we chose as the 
pivot scale k = 0.002 Mpc" 1 (Peiris et al. 2003). 

We perform a parameter analysis as described in 
Section 3.1 at this pivot scale with the addition of ten- 
sor modes with a free amplitude ratio r < 1 and the 
constraint in equation (11). In Figure 4 we show the re- 
sults of this analysis. The plot shows the results in the 
Vv-^v plane. The models are colour coded according to 
the third slow-roll parameter ev- From equation (9) we 
see that < ev < 0.0625, if we restrict the analysis 
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to r < 1. For most models of inflation there is a negli- 
gible contribution of gravitational waves to the pertur- 
bations, i.e. r ~ ev « (Liddle & Lyth 2000). In this 
case the likelihood contour in Figure 4 is restricted to 
the area depicted by the dark points, towards the lower 
left of the figure, r/v and £v are thus relatively small, 
which according to equations (8) and (10) corresponds 
to a spectrum very close to scale invariant scalar fluc- 
tuations. This is consistent with our findings in Section 
3.1. However, with the current data, models with a large 
tensor contribution are also valid, which correspond to 
the light coloured points in Figure 4, in the top right 
corner. The marginalized best fit values for the slow-roll 
parameters are e v = 0.032 ± 0.018, r/ v = 0.020 ± 0.025 
and £v = 0.015 ±0.014. The constraint on the tensor to 
scalar ratio is r = 0.499 ± 0.296. 



4 CONCLUSION 

We have combined the linear dark matter power spec- 
trum inferred from the Lyman-a forest with the 1st year 
WMAP data to obtain constraints on the rms fluctu- 
ation amplitude of matter fluctuations as and on the 
slope and shape parameter of the primordial density 
power spectrum n s and n run . We have restricted our 
analysis to spatially flat models with a low gravita- 
tional wave contribution to the CMB fluctuations. The 
WMAP data alone provides only weak constraints on 
erg, n s and n run . Adding the Lyman-a forest significantly 
tightens the constraints because of the resulting long 
lever arm which spans more than 3 orders of magni- 
tude in wavenumber. The WMAP data combined with 
the DM power spectrum inferred by Croft et al. (2002) 
from their Lyman-a forest data requires a tilt (n s < 1) 
and/or running of the spectral index (n run < 0) to fit 
the data which is significant at the 3 a level. If we com- 
bine instead the WMAP data with the DM power spec- 
trum inferred by VHS a smaller tilt and running of the 
spectral index are required which are not statistically 
significant. The best fit value for the rms fluctuation 
amplitude of matter fluctuations is cr 8 = 0.94 ±0.08 and 
n s = 0.99 ± 0.03, in good agreement with the simpler 
analysis presented by VHS. 

The difference of the analysis for the DM power 
spectrum of C02 and VHS is mainly due to the more 
appropriate smaller effective optical depth obtained 
from high resolution, high S/N QSO absorption spec- 
tra which was assumed in VHS. We further included 
tensor modes and put constrains on the slow roll param- 
eters. At present the limits on ev, (,v and r\y are not 
very tight, but future CMB missions such as Planck, 
together with better understanding of the systematic 
uncertainties of the Lyman-a forest data will improve 
the constraints. 
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